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Abstract. In this paper total cross sections for signals and backgrounds of top- and Higgs-production 
channels in e^e~ collisions at a future linear collider are presented. All channels considered are characterized 
by the emergence of six-particle final states. The calculation takes into account the full set of tree-level 
amplitudes in each process. Two multi-purpose parton level generators, HELAC/PHEGAS |1I2| and AMEGIC++ 
12] are used, their results are found to be in perfect agreement. 



1 Introduction 

Six-particle final states constitute the signature for many 
processes that will be studied at the precision level at a fu- 
ture e"'"e~ linear collider. Important channels include the 
production and decay of top quark pairs and - if existent - 
of one or more Higgs bosons, the latter process allowing a 
test of the structure of the Higgs potential. Furthermore, 
if no evidence for a Higgs boson was found at the LHC, 
the study of quartic gauge boson couplings is mandatory 
in order to understand alternative scenarios of electroweak 
symmetry breaking. Leaving the framework of the Stan- 
dard Model (SM) the production of, say, chargino pairs 
in the Minimal Supersymmetric Standard Model (MSSM) 
will lead to six-particle final states as well. To understand 
these processes at the precision level, i.e. at the order of 
a few per cent, it is mandatory to supplement typical ap- 
proaches such as the narrow-width approximation, with 
corresponding calculations through full amplitudes, and 
to quantify the effect of non-resonant contributions. Ob- 
viously, for hadronic final states, a full QCD calculation 
is unavoidable. 

Such investigations, however, are a formidable calcula- 
tional task that cannot be handled without dedicated com- 
puter programs. Two major difficulties make these neces- 
sary: 

1. Including the full SM for the production of a six-particle 
final state often leads to having to handle a large num- 
bers of diagrams. As an illustrative example of this 
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problem, take the process e+e^ e^e"e+e^e+e^, 
which results in 13896 Feynman diagrams. Obviously, 
the common textbook method of squaring the dia- 
grams, by employing completeness relations for the ex- 
ternal particles and evaluating the traces, is not a very 
efficient way to calculate the matrix element squared. 
2. Apart from the treatment of an enormous number of 
diagrams, growing roughly factorially with the num- 
ber of external legs, the integration over the phase 
space of the outgoing particles becomes a tedious task. 
The high dimensionality, 3n — 4, for n final-state par- 
ticles necessitates the use of Monte Carlo methods. 
To achieve convergence of the Monte Carlo procedure 
process- and cut-dependent phase-space mappings are 
required that tame wildly fluctuating integrands, which 
are due to nearly on-shell propagators. A benefit of 
Monte Carlo methods, if carefully implemented, is that 
not only total cross sections but also distributions in- 
cluding all possible types of kinematical cuts can be 
calculated on an equal footing. 

In the past years, different types of parton level generators 
have been constructed. They can be crudely characterized 
as either specialized or multi-purpose generators. 
Usually, the former ones contain explicit matrix elements 
and phase-space mappings for specific classes of processes 
with specific assumptions. These matrix elements were 
constructed before outside the respective program, and 
this feature also allows for instance to implement non- 
universal higher order corrections in a controlled way. Of- 
ten, the phase space mappings can be optimized before as 
well. 

Examples for such programs dealing with some of the pro- 
cesses discussed here are LUSIFER 0] and eett6f jS]. Both 
are constrained to fermions in the final state; in the case of 
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LUSIFER these fermions are bound to be massless, whereas 
eett6f specializes in top quark pair production channels 
where the outgoing fermions might be massive, but elec- 
trons are disallowed in the final state. Both programs use 
versions of the adaptive multichannel method in the spirit 
of [H| for their integration. A further dedicated program 
using the multichannel importance sampling is SIXFAP 7 . 
It provides the electroweak contributions for a large set of 
six-fermion final states, taking into account possible non- 
zero fermion masses. 

In contrast to specialized programs, multipurpose codes 
generate both the matrix elements and the phase space 
mappings with or without some intervention by the user. 
Apart from the programs used in this paper, examples 
of these types of programs are D 'Mega/Whizard 8 ,9 and 
MadGraph/MadEvent [TUIITT] . In the first package, O'Mega 
|H] relies on a version of the alpha algorithm How- 
ever, in the present version of O'Mega, full QCD has not 
been implemented yet. The integration of the resulting 
matrix elements is achieved through Whizard jQj, which 
constructs phase-space mappings automatically and inte- 
grates them with the VAMP-algorithm In fact, Whizard 
can also be interfaced with other matrix element gen- 
erators and it can be used to generate unweighted, sin- 
gle events. In contrast, MadGraph/MadEvent generates all 
Feynman diagrams for a process under consideration and 
then passes the information to the HELAS package [T^ 
for the translation into corresponding helicity amplitudes. 
The integration channels are constructed automatically, 
and a new version of the adaptive multichannel method 
described in JTj is employed for the actual integration and 
the generation of unweighted events. 

This paper deals with the calculation of total cross sec- 
tions for many relevant processes at a future e'^e~ col- 
lider with two different, independent packages, namely 
HELAC/PHEGAS and AMEGIC++. Similar to the comparison 
of four-fermion generators at the LEP2 Monte Carlo Work- 
shop |15| . a detailed study and mutual benchmarking of 
tools for six- and eight-particle final states at a future 
linear collider has been initiated in the framework of the 
extended ECFA/DESY study [El. Here, a further step 
into this direction is reported. 

For the case of only massless final state particles, a 
similar comparison between the programs LUSIFER and 
MADGRAPH, the latter using WHIZARD for the phase-space 
integration has been presented in 4 . In addition, results 
achieved by different generators for selected top quark pair 
production channels in the massless fermion approxima- 
tion can be found in |T7] . 

The paper is organized as follows. In Sect. [5] and |21the 
relevant features of the two programs, HELAC/PHEGAS and 
AMEGIC++, are briefly reviewed. In Sect.0]the results are 
presented and discussed. Conclusions are drawn in Sect.^l 



2 The HELAC/PHEGAS package 

2.1 Amplitude Computation: HELAC 

The traditional representation of the scattering amplitude 
in terms of Feynman graphs results in a computational 
cost that grows like the number of those graphs, therefore 
as 7i!, where n is the number of particles involved in the 
scattering process. 

An alternative to the Feynman graph representation is 
provided by the Dyson-Schwinger approach. Dyson-Schwinger 
equations recursively express the n-point Green's func- 
tions in terms of the l-,2-,. . . , (n — l)-point functions. For 
instance in QED these equations can be written as follows: 




n 

1=1 

+ {^aW^{P2)i''i^{PMPi,P2) (2) 

P=Pl+P2 

where 

h^{P) = ^{P) = ^ ^(P) = ^ (3) 

describes a generic n-point Green function with one out- 
going photon, fermion or antifermion leg, respectively, car- 
rying momentum P. Flf^u stands for the boson propagator 
and e takes into account the sign due to fermion anti- 
symmetrization. 

Technical details about the implementation of the algo- 
rithm for the electroweak interactions can be found in . 
For QCD amplitudes, colour representation and summa- 
tion play an important role. Usually, for the n-gluon am- 
plitude the well known colour decomposition is used 

M = 2ig''-^ Tr(i''i ...t'^")C(l,...,n) (4) 

P(2,...,n) 

For processes involving quarks a similar expression may be 
derived. For further details, the reader is referred to the 
vast literature on the subject ^H]- Methods for calculating 
the C-functions have been developed JH| j including some 
recent ones, more suitable for multiparticle processes [201 
I21| . One of the most interesting aspects of this decomposi- 
tion is the fact that the C-functions satisfy certain useful 
properties, such as gauge invariance and cyclic symmetry. 
Nevertheless, the computational complexity is rather high 
and the evaluation of the squared colour matrix a rather 
complicated task |22| . 

In HELAC a novel approach has been considered. It is based 
on the colour connection (or colour flow) representation of 
the interaction vertices, where the explicit reference to the 
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colour has been avoided, as is also the case in the usual 
colour decomposition. The advantage, however, is that the 
colour factors aquire a much simpler form, which moreover 
holds for gluon as well as for quark amplitudes, leading 
to a unified approach for any tree-order process involving 
any number of coloured partons. Moreover, the unweight- 
ing procedure is significantly facilitated, since the usual 
information on colour connections, needed by the parton 
shower Monte Carlo, is automatically available, without 
any further calculation. The colour factor is universally 
given by 

= ^lai(l)S2<j,i2) ■ ■ ■ Sncr,{n) j (5) 

whereas the colour matrix, defined as 

colours 

with the summation running over all colours, 1 . . . Nc, has 
a very simple representation: 

Mij = Ar™('^^^'^-') . (7) 

Here, 1 < m(ai,<7j) < n counts the number of cycles 
made by the elements of the permutations tr/ and aj. 
Details can be found in ref. |23| . 

2.2 Phase-space integration: PHEGAS 

The study of multiparticle processes, such as six-fermion 
production in e^e^, requires efficient phase-space Monte 
Carlo generators. The reason is that the squared ampli- 
tude, being a complicated function of the kinematical vari- 
ables, exhibits strong variations in specific regions and/or 
directions of the phase space, lowering in a substantial 
way the speed and the efficiency of the Monte Carlo inte- 
gration. A well known way out of this problem relies on 
algorithms characterized by two main ingredients: 

1. The construction of appropriate mappings of the phase 
space parametrization, in such a way that the main 
variation of the integrand can be described by a set of 
almost uncorrelated variables, and 

2. A self-adaptation procedure that reshapes the gener- 
ated phase-space density in order to be as close as pos- 
sible to the integrand. 

In order to construct appropriate mappings, it is impor- 
tant to note that the integrand, i.e. the squared amplitude, 
has a well-defined representation in terms of Feynman di- 
agrams. It is therefore natural to associate to each Feyn- 
man diagram a phase-space mapping that parametrizes 
the leading variation coming from it. In PHEGAS, infor- 
mation from HELAC is used to automatically construct a 
representation of all Feynman graphs contributing to the 
given process. The subset of Feynman graphs that results 
in a different phase-space parametrization is then used 
as kinematical mappings, called channels, to perform the 
Monte Carlo integration. Details can be found in |2j. 



Since in six- and eight-fermion production a large number 
of kinematical channels contribute, typically of the order 
of 10^ to 10"*, the optimization is also used to reduce their 
number. This is based on the fact that many of the chan- 
nels exhibit an important correlation that renders them 
practically useless as separate channels. The reduction in 
the number of channels achieved by this optimization is 
generally important, resulting in a very efhcient and rapid 
integration. 

The main points can be summarized as follows: 

— The algorithm exhibits a computational cost that grows 
like ^ 3", in contrast to the nl growth of the Feynman 
graph approach. Therefore there is no severe limitation 
in computing many-particle amplitudes (up to at least 
12 external). 

— All electroweak vertices in both the Feynman and the 
Unitary gauge have been included, allowing highly non- 
trivial checks to be performed. The QCD interactions 
have been implemented in the colour-connection repre- 
sentation, allowing also a fast unweighting procedure. 
The decay width of unstable particles is introduced in 
the fixed- width and complex-mass schemes. Any pro- 
cess with any type of Standard Model particle can be 
reliably computed. 

— Special features include also the possibility to use higher 
precision floating point arithmetic, allowing full con- 
trol over all possible phase-space regions. Speeding up 
techniques, for hclicity Monte Carlo treatment and 
large Nc estimates, are also available. 

— Incorporation of higher order corrections (currently 
available Fermion-Loop corrections up to three-point 
functions) and the introduction of the Minimal Super- 
symmetric Standard Model are in progress. 

3 The program AMEGIC++ 

AMEGIC++, acronym for (A Matrix Element Generator in 
CH — h), is a multi-purpose parton- level generator written 
in C++. It provides a convenient tool for the calculation 
of cross sections for scattering processes at the tree level 
in the framework of the SM and the MSSM. Recently the 
code was extended to cover processes in the ADD model 
of large extra dimensions as well '24 . The program can 
also be used to generate single events and it is one of the 
modules for the new complete event simulation framework 
SHERPA As such, the single events of AMEGIC++ can be 
handed over to the parton shower module APACIC++ 
with the help of a new method that is correct at the next- 
to-leading logarithmic accuracy [211 ^^^^ linked 
correctly to fragmentation. 

In AMEGIC++, full sets of Feynman diagrams are constructed 
automatically and are translated by the program into he- 
licity amplitudes in a formalism similar to the one in |28| . 
The colour structure of each diagram is represented as 
a word string, the emerging structures are grouped into 
sets of amplitudes with identical, common colour struc- 
ture. Based on them, a matrix of colour factors between 
amplitudes is calculated using the ordinary SU{3) alge- 
bra. A number of refinements of the helicity method has 
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been implemented within the code as well. First of all, 
the algorithm presented in fixes the relative signs 
of amplitudes when Majorana fermions are present. Fur- 
thermore, explicit polarizations for massive or massless 
external spin-1 bosons are enabled, allowing us to con- 
sider polarized cross sections. Similar considerations help 
to replace numerators of spin-1 propagators by summing 
over suitably defined polarizations for off-shell particles 
disentangling nested Lorentz structures emerging for am- 
plitudes with many internal spin-1 bosons. As a result, 
AMEGIC++ needs only quite a limited set of building blocks 
to construct all helicity amplitudes. Internally, they are 
represented as word strings employing some knowledge- 
storing mechanism that ensures that all building blocks 
have to be evaluated only once for each call of the full 
matrix element. With the help of internal methods these 
word strings are further simplified. Furthermore, another 
massive gain in efficiency has been achieved by summing 
amplitudes with identical colour structure and by algo- 
rithms for finding common factors. This is exemplified in 
Fig.n Having performed these manipulations, the result- 
ing helicity amplitudes are stored in library files. 




Fig. 1. Factoring out common pieces of amplitudes with iden- 
tical colour structure. In the example above, the parts within 
the boxes are identical, hence the two amplitudes can be added 
and the terms inside the box can be factored out. 

There are a number of prescriptions to treat unstable par- 
ticles. At the moment, AMEGIC++ supports the fixed-width 
scheme (FWS) and the complex-mass scheme (CMS). 
Defining the complex mass parameters of the electroweak 
gauge bosons, the Higgs boson and the top quark in terms 
of the real masses and the constant widths through 

^1 = mjf - iFh niH , Mt = mt - iTtll , (8) 
the corresponding propagators can be written as 

In the FWS, the electroweak mixing angle is defined ac- 
cording to 

sin2 ew = l-'^. (10) 

It is kept real. For the case of the gauge-invariant CMS, 
the real gauge-boson masses have to be replaced by their 



complex counterparts and this parameter is therefore com- 
plex as well. 

Within AMEGIC++ the Yukawa couplings of fermions to the 
Higgs boson and their kinematical masses are decoupled. 
This allows us to study, for example, the production of 
Higgs bosons and their decay into 5-quarks, even in those 
cases where the user prefers to neglect the influence of the 
6-mass on both the phase space and the helicity structure. 
For the integration over the phase space of the outgo- 
ing particles, AMEGIC++ employs an adaptive multichan- 
nel method 6 . Similar to their implementation, generic 
elements for phase-space mappings such as propagator- 
like structures are provided. The individual Feynman dia- 
grams are analyzed individually and one or more suitable 
phase-space parametrizations for each diagram are auto- 
matically created and stored in library files. As an ex- 
ample, consider Fig. EJ which exhibits a diagram and its 
translation into propagator- and decay-parametrizations. 
These files, both for the amplitudes and the phase-space 
parametrizations, are compiled and linked to the code be- 
fore the actual integration starts. 
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Fig. 2. Translation of a Feynman diagram into a phase-space 
parametrization. Ds,a denote symmetric or asymmetric decays 
; the latter ones reproduce the typical feature of coUinear emis- 
sion of particles notorious for gauge theories with massless 
spin-1 bosons. The propagator terms for massless particles Pq 
peak at the minimal allowed invariant mass. 

For users of AMEGIC++ only very little intervention is needed. 
Having specified the process(es), the model framework and 
its parameters, a first "initialization" run of the code re- 
sults in the creation of library files. After their compila- 
tion, a second, "production" run will generate the results 
without any further manipulation. 



4 Numerical results 

4.1 Input parameters and phase-space cuts 

The SM parameters are given in the scheme: 

mw = 80.419 GeV , Fw = 2.12 GeV, 
mz ^ 91.1882 GeV , Fz = 2.4952 GeV, 
Gf, = 1.16639 X 10"^ GeV"^ 
sin^ew = 1 - mly/m%, 

as 0.0925(0.0891) at 360(500) GeV. (11) 
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The electromagnetic coupling is derived from the Fermi 
constant according to 

TT 

The mass of the Higgs boson is assumed to he Mh = 
130 GeV and its associated SM tree-level width is Fh — 
0.00429 GeV. For this Higgs boson mass its branching ra- 
tios H ^ bb and H — > W'^ W~ 4/ are of the same 
order and therefore both decay channels signify its occur- 
rence as an intermediate state. For the massive fermions, 
the following masses have been used: 

= 105.6583 MeV, = 1.777 GeV, 
TO„ = 5 MeV, nid = 10 MeV, 
TO, = 200 MeV, TOc = 1.3 GeV, 
TOh = 4.8 GeV, 

TOt = 174.3 GeV, Ft ^ 1.6 GeV. (13) 

The constant widths of the electroweak gauge bosons, the 
Higgs boson and the top quark are introduced via the 
fixed- width scheme as defined in Sect. 13 CKM mixing of 
the quark generations and the coupling of the Higgs boson 
to the very light fermion flavors (e, u, d) is neglected. 
Concerning the phase-space integration, the following cuts 
are applied on the external particles: 

0{l, beam) > 5° , 0{l, I') > 5° , > 10 GeV, 
d{q, beam) > 5° , 0{l, (?) > 5° , > 10 GeV, 
m{q,q') > 10 GeV , (14) 

where 0{i,j) specifies the angle between the particles i and 
j in the center-of-mass frame, and I, q and beam denote 
charged leptons, quarks or gluons and the beam electrons 
or positrons, respectively. The invariant mass of a jet pair 
qq' is denoted by m{q,q'). 

All results presented here are obtained using 10^ points 
(before cuts); statistical errors of the Monte Carlo inte- 
grations, i.e. one standard deviation, are given in paren- 
theses. 



4.2 Results 

First of all, processes have been considered that serve as 
signals or backgrounds for the production and decay of 
top pairs. Tabled Since the branching ratio is practically 
100% for the decay of_top quarks into bottom quarks and 
a.W [t ^ bW^, i — > bW~), all modes considered include 
a pair of bottom quarks. In cases involving a mixture of 
top production and decay and pure QCD diagrams, the 
relative importance of the different contributions to the 
total cross section has been estimated by switching on 
and off the QCD coupling constant. In both cases (the 
fully hadronic mode bbuudd and the semileptonic mode 
bbude~De) the top contribution is by far the dominating 
channel; the difference of taking into account the QCD 
contributions or neglecting them is of the order of 2-3%. 



Top-quark channels 



Final state 


QCD 


AMEGIC++ [fb] 


HELAC [fb] 


- - 
bbuddu 


yes 


32.90(15) 


33.05(14) 




yes 


49.74(21) 


50.20(13) 




no 


62.ZZ[64: ) 


oo 1 o / 1 n\ 
32.12(19) 




no 


A r\ A / A A \ 

49.42(44) 


50.55(26) 


bouugg 




11.23(10) 


11.1od(41) 






n 1 1 / 1 o^ 

y.ii(ioj 


(i.o6z[4:6) 


- 

bbgggg 




10 oo/ioN 
18.82(13) 


1 O TA / 1 1 \ 

18.79(11) 






24.U9(iB) 


oo on/' 1 '7\ 
23.80(1 ( ) 


bbude~Ve 


yes 


11.460(36) 


11.488(15) 




yes 


17.486(66) 


17.492(41) 




no 


11.312(37) 


11.394(18) 




no 


17.366(68) 


17.353(31) 


bbe^ Ve^ Ve 




3.902(31) 


3.885(7) 






5.954(55) 


5.963(11) 


bbe'^UeiJ.^i'f^ 




3.847(15) 


3.848(7) 






5.865(24) 


5.868(10) 


bbn'^i/^jj.^i'^ 




3.808(16) 


3.861(19) 






5.840(30) 


5.839(12) 



Table 1. The cross sections for possible signals and back- 
grounds of top quark pair production in e~ annihilation. 
All results in fb for = 360 GeV (first row) and = 500 
GeV (second row). 



Also, the total cross section of the fully hadronic channel 
is substantially larger than the cross section of any other 
individual bb+A jets mode. 

For the QCD contributions, a similar pattern arises 
also in the vector-boson fusion channels, cf. Tables El and 
121 These channels are characterized by either an electron- 
positron or an electron-neutrino anti-neutrino pair in the 
final state, corresponding to either Z boson or to W boson 
fusion processes, respectively. Again, switching on and off 
the QCD coupling constant gives rise to differences on the 
level of a few per cent. In contrast, taking into account the 
Higgs boson (Table O which may be produced in the s- 
channel through the fusion of two t-channel vector bosons, 
or neglecting it (Table O changes the total cross sections 
for all channels considered by a factor of 2 or larger. This 
is especially pronounced for channels that can be identi- 
fied as VFW^-fusion channels with a semileptonic or fully 
hadronic decay of the W^-pair_produced by the Higgs decay 
(i.e. Vf.i>eude~i'e and Vf.i>eudii~v^, or UeVeuddu, respec- 
tively) , where the cross sections are larger by one order of 
magnitude. 

Another mode for Higgs production at an electron- 
positron colhder is Higgs-strahlung, where the Higgs bo- 
son is radiated off a Z-boson in the s-channel. In TableQ] 
total cross sections for such modes are displayed, where 
the Z boson decays into muons and the Higgs boson goes 
into four fermions through a pair of W ot Z bosons. In 
Table |S1 identical total cross sections for the same final 
states, but neglecting the Higgs contribution, are shown. 
In both cases, again, the size of the pure QCD contri- 
butions is found to be small for most final states, i.e. of 
the order of few per cent. The only exception is for a 
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Vector fusion with Higgs exchange Higgs production through Higgs-strahlung 



Final state 


QCD 


AMEGIC++ [fb] 


HELAC [fb] 


e e uudd 


yes 


U. 0042(85) 


U.d858(o1) 




yes 


l.z37(15) 


1.265(5) 




no 




U.OOZ/(dO) 




no 




i.Zo94( / 0) 


e e uue e 




o.\)o(Sb)c-OS 


D.113(57)e-03 






D.58(23)e-U3 


a G'\ A ( Of\ \ ^ AO 

D.ol4(50)e-03 


- + - - + 
e e uuji n 




y.z4(lz je-Uo 


n A/l/'1^^.^ AQ 
y.U4(ll )e-Uo 






9.25(17)e-03 


9.145(74)e-03 


VeVeUddu 


yes 


1.15(3) 


1.176(6) 




yes 


2.36(7) 


2.432(12) 




no 


1.14(3) 


1.134(5) 




no 


2.35(7) 


2.429(13) 


Vei'eUde' i>e 




0.426(11) 


0.4309(48) 






0.916(30) 


0.9121(48) 


VeVeUdjl'Dij, 




0.425(12) 


0.4221(30) 






0.878(27) 


0.8888(47) 



Table 2. The cross sections for different e e~ ^ 6/ final 
states corresponding to the Higgs production via vector-boson 
fusion signal. All results in fb for ^/s = 360 GeV (first row) 
and = 500 GcV (second row). 

Vector fusion without Higgs exchange 



Final state 


QCD 


AMEGIC++ [fb] 


HELAC [fb] 


e~e^uudd 


yes 


0.4838(50) 


0.4842(25) 




yes 


1.0514(97) 


1.0445(51) 




no 


0.4502(31) 


0.4524(23) 




no 


1.0239(79) 


1.0227(43) 






3.757(98)e-03 


3.577(43)e-03 






4.082(56)6-03 


4.214(46)e-03 






5.201(61)e-03 


5.119(70)e-03 






5.805(67)e-03 


5.828(49)e-03 


Ve^eUddu 


yes 


0.15007(53) 


0.15070(64) 




yes 


0.4755(21) 


0.4711(24) 




no 


0.12828(42) 


0.12793(55) 




no 


0.4417(19) 


0.4398(21) 


VeVeUde' Ve 




0.04546(13) 


0.04564(19) 






0.16033(63) 


0.16011(78) 


VeVeUdjl'Vij, 




0.04230(12) 


0.04180(16) 






0.14383(53) 


0.14439(65) 



Table 3. The backgrounds to Higgs production via vector 
boson fusion. All contributions from intermediate Higgs bosons 
are neglected. Cross sections are given in fb for ^/s = 360 GeV 
(first row) and ^/s = 500 GeV (second row). 



pair of muons and four identical quarks; there, the inclu- 
sion of QCD changes the results by roughly 20%, when 
the Higgs boson is taken into account, and by a factor of 
roughly 2 when its contribution is neglected. It is amusing 
to note that this relative factor of two compares in size 
with the effect of including the Higgs boson itself. This, 
however, is true only for the mode that can be imagined as 
e^e~ — > ZH — > ZZZ ii^ fi^ uuuu. In all other cases, as 
said before, inclusion of QCD has minor effects only; the 



Final state QCD AMEGIC++ [fb] HELAC [fb] 





>^ 




- 


0.03244(27) 


0.03210(15) 








- 


0.03747(29) 


0.03749(32) 


At" 


>^ 


^ude~i>e 


- 


0.0924(8) 


0.09306(46) 










A 1 1 A^? /'00\ 

0.1106(22) 


A "1 AAAT f c:c:\ 

0.10901(66) 




>^ 


^ fj,~ fj,'^ e~ e'^ 




2.828(67)c-03 


2.923(52)c-03 










2.731(65)e-03 


2.691(42)6-03 




>^ 


^uudd 


yes 


0.2534(24) 


0.2540(16) 








yes 


0.2634(22) 


0.2642(15) 








no 


0.2441(23) 


0.2471(15) 








no 


0.2593(22) 


0.2589(14) 




>^ 


'uuuu 


yes 


1.125(8)e-02 


1.135(22)c-02 








yes 


8.767(65)e-03 


8.978(58)e-03 








no 


7.929(57)e-03 


8.078(92)c-03 








no 


6.098(35)6-03 


6.013(26)c-03 



Table 4. The cross sections for different e e~ — > 6/ final 
states corresponding to the Higgs-strahlung signal. All results 
given in fb for •Js = 360 GeV (first row) and = 500 GeV 
(second row). 

Backgrounds to Higgs-strahlung 



Final state 


QCD 


AMEGIC++ [fb] 


HELAC [fb] 










0.01845(14) 


0.01843(13) 










0.03054(23) 


0.03092(19) 




-At^ 


^ude~Ve 




0.05284(57) 


0.05209(33) 










0.08911(53) 


0.08925(48) 






^ fj,~ fj,'^ e~ e'^ 




2.204(52)c-03 


2.346(49)c-03 










2.280(66)c-03 


2.277(62)c-03 


At" 


-At^ 


^uudd 


yes 


0.1412(10) 


0.1404(11) 








yes 


0.2092(12) 


0.2075(13) 








no 


0.1358(20) 


0.1341(12) 








no 


0.2040(12) 


0.2015(11) 


A^" 




^uuuu 


yes 


5.937(24)c-03 


5.937(25)c-03 








yes 


6.134(29)6-03 


6.108(27)c-03 








no 


2.722(10)6-03 


2.710(11)6-03 








no 


3.290(12)6-03 


3.303(12)6-03 



Table 5. Background contributions to the Higgsstrahlungs 
signal for various 6/ final states. All diagrams with intermedi- 
ate Higgs bosons have been neglected. Cross sections are given 
in fb for = 360 GeV (first row) and = 500 GeV (second 
row). 



Higgs boson in contrast roughly doubles the total cross 
section in all the other channels. 

One of the salient research goals at a potential linear 
collider operating at energies around 500 GeV is the deter- 
mination of the Higgs potential. For this, the self- couplings 
of the Higgs bosons have to be checked. In the framework 
of this publication, results arc provided for the channel 
where the Higgs bosons emerge in Higgs-strahlungs-like 
topologies and decay into a pair of bottom quarks. This 
leads to final states + 4&, where the muons mainly 

come from the Z bosons. Results for total cross sections for 
the process e+e" IJ^^ + 4&, where contributions me- 
diated by Higgs bosons have been included or neglected, 
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Triple 


Higgs coupling 




Final state QCD 


AMEGIC++ [fb] 


HELAC [fb] 


jj.^ jj.'^bbbb yes 


2.560(26)e-02 


2.583(26)e-02 


yes 


3.096(60)6-02 


3.019(43)6-02 


no 


1.711(55)e-02 


1.666(28)6-02 


no 


2.34(12)e-02 


2.36(10)6-02 



Table 6. Cross sections for the process e"*" e ^ A* iJ.'^bbbb. 
All results in fb for ^ = 360 GeV (first row) and ^/s = 500 
GeV (second row). 

Backgrounds to triple Higgs coupling 
Final state QCD AMEGIC++ [fb] HELAC [fb] 



fi fi'^bbbb 


yes 


7.002(32)6-03 


7.044(22)6-03 




yes 


6.308(24)6-03 


6.364(21)6-03 




no 


2.955(11)6-03 


2.972(12)6-03 




no 


3.704(15)6-03 


3.695(13)6-03 



Table 7. Cross sections for e ~* j-i jj.'^bbbb with all contri- 
butions due to intermediate Higgs bosons left out. All results 
in fb taken for ^/s = 360 GeV (first row) and = 500 GeV 
(second row). 



are given in Tables and |7| respectively. From the results 
displayed one can read off that the inclusion of intermedi- 
ate Higgs bosons enhances the cross sections by a factor 
of three to four. Again, also the effect of QCD has been 
checked. For the process involving the intermediate Higgs 
bosons, QCD leads to total cross sections that are larger 
by roughly 30%-40%, without the Higgs bosons, QCD con- 
tributes on the level of factors of two to three. 



5 Summary of results 

In the framework of this comparison, total cross sections 
for 86 different processes involving six-particle final states 
have been obtained by the two multi-purpose matrix el- 
ement generator packages HELAC/PHEGAS and AMEGIC++. 
The integration over the multidimensional phase space 
of the final states has been performed with Monte Carlo 
methods, and in all cases one million MC points have been 
used. For nearly all cross sections the resulting statistical 
error was significantly smaller than one per cent, roughly 
five per mille. There have been no significant differences 
between the two codes. To compare the results, for each 
process i the deviation s^'^ of the two resulting cross sec- 
tions tJy'' and cr^'' has been calculated through 



(15) 



Aa 



ii) 



Aa 



The distribution of the individual differences is depicted in 
Fig. El The average deviation is s = —0.065, the variance 
in their distribution is ctc w 1. The maximal difference 



HELAC/PHEGAS vs. AMEGIC++ 

differences in results 

' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 



average ; -0.065, variance : 1 



10 - 



IT 



statistical deviation in o 

Fig. 3. The distribution of deviations s'*', given by Eg 1151 for 
the eighty-six total cross sections i presented in this paper. The 
average value is s = —0.065, their variance is as ~ 1 . 



between two cross sections is smaller than three standard 
deviations, gC™'*'^ ) 2.6. The distribution of differences 
follows roughly a Gaussian distribution. 

To summarize: Both packages, HELAC/PHEGAS as well 
as AMEGIC++, lead, with quite different methods, to con- 
sistent results for total cross sections for a large number 
of different processes with six particles in the final state. 
This provides an independent check of the precision level 
of the two codes, which can be considered as successfully 
tested. 
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